NPS-MA-99-001 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


CONVOLUTION  METHODS  FOR 
MATHEMATICAL  PROBLEMS  IN  BIOMETRICS 

by 

Christopher  Frenzen 
January  15,  1999 


Approved  for  public  release;  distribution  is  unlimited. 

Prepared  for:  Biometrics  Identification 
College  of  Engineering 
San  Jose  State  University 
One  Washington  Square 
San  Jose,  CA  95192-0205 


1  99902180  84 


REPORT  DOCUMENTATION  PAGE 

. 

Form  approved 

0MB  No  0704-0188  . 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources. 
Gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  information  Operations  and  Reports,  1215  Jefferson  Davis 
Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  OfBce  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-01S8),  Washington,  DC  20503. 

1.  AGENCY  USE  ONLY  (Leave  blank) 

2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

January  15, 1999  March  1997  -  September  1997 

4.  TITLE  AND  SUBTITLE 

Convolution  Methods  for  Mathematical  Problems  In  Biometrics 

5.  FUNDING 

MIPR  #  NAFRLP62800493 

6.  AUTHOR(S) 

Christopher  L.  Frenzen 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School 

Monterey,  CA  93943 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

NPS-MA-99-001 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Biometrics  Identification 

College  of  Engineering 

San  Jose  State  University 

One  Washington  Square 

San  Jose,  CA  95192-0205 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

Approved  for  public  release;  distribution  is  unlimited. 

12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words.) 

Estimation  of  the  impostor  probability  density  function  from  the  inter-template  and  template-sample 
probability  density  functions  involves  a  multi-dimensional  convolution  of  the  latter  density  functions.  We 
assume  isotropy  of  the  probability  distribution  functions  and  use  the  Fourier  transform  to  express  the 
convolution  as  a  one-dimensional  integral  with  a  kernel  given  by  a  spherical  Bessel  function. 

14.  SUBJECT  TERMS 

Biometrics,  Convolution 

15.  NUMBER  OF 

PAGES 

6 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION 

OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

UNCLASSIFIED  UNCLASSIFIED  UNCLASSIFIED 

20.  LIMITATION  OF 
ABSTRACT 

NSN  7540-01  -280-5800  Standard  Form  298  (Rev.  2-89) 


Prescribed  by  ANSI  Std  239-18 


Enclosure  (6) 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California  93943-5000 


RADM  Robert  C.  Chaplin  R.  Elster 

Superintendent  Provost 

This  report  was  prepared  for  Biometrics  Identification,  College  of  Engineering,  San  Jose  State 
University. 

This  report  was  prepared  by: 


Christopher  L.  Frenzen 
Professor  of  Mathematics 


Reviewed  by: 


Department  of  Mathematics 


Released  by: 


Dean  of  Research 


Convolution  Methods  for  Mathematical  Problems  in  Biometrics 


C.L.  Prenzen 

Department  of  Mathematics 
Naval  Postgraduate  School,  Monterey,  CA 


1. Introduction.  The  problem  we  shall  investigate  can  be  formulated 
in  the  following  way,  due  to  Peter  Bickel  [1]:  For  M  +  1  individuals  K 
measurements  Xij,  i  =  1,  •  •  • ,  M  +  1,  j  =  1,  •  -  • ,  iiT  are  made.  The  Xij  are 
assumed  to  be  vectors  in 

We  assume,  following  Bickel  [1],  that 

Xij  fXi  "h  Sij, 

where  /^i,  •  •  • ,  f^M+i  are  individual  effects  vectors  from  which  are  i.i.d.  (in¬ 
dependently  and  identically  distributed)  with  a  distribution  which  is  spheri¬ 
cally  symmetric  about  some  jx  in  R^.  The  are  i.i.d.  independently  of  the 
Hi  and  their  distribution  is  spherically  symmetric  about  0. 

The  K{M  -h  l)A'ijS  represent  K  biometric  measurements  made  on  M  -fl 
individuals;  however  the  measuring  device  actually  records  only  the  data 
I  Xi  —  Xi!  I  and  |  A'm+ij  —  Xm-^-i  |,  1  <  i  <  Af,  1  <  j  <  fc,  where  |  x  |  is 
the  length  of  x  and 

1  ^ 
y-i 

is  the  centroid  of  the  K  measurements  taken  on  the  the  zth  individual. 

Let  p  be  the  density  oi  \  Xi  —  X2  \,  q  the  density  of  |  Xm+x,\  —  Xm+i  |, 
and  r  the  density  of  |  Xm+i,i  I-  The  basic  problem  is  to  estimate  r  given 
estimates  p,  q  for  p  and  q.  As  p  and  q  are  densities  on  [0, 00),  the  problem  of 
estimating  p  and  g  is  a  standard  problem  in  estimation  theory  which  we  do 
not  consider  further. 


2.  Bickel’s  Approach.  What  is  the  relationship  between  p,  q,  and  r? 
Following  Bickel  [1],  we  write 

Xm+1,1  —  X\  =  Sm+i,!  +  hm+\  “  Ml  —  £1)  (1) 
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~  ^M+1  +  fJ-M+l  ~  (J'l—  (2) 

^M+l,l  —  ^M+1  =  ^M+1,1  -~  £m+1-  (3) 

Note  that  addition  of  equations  (2)  and  (3)  yields  equation  (1).  Further, 
the  quantities  on  the  left  sides  of  (2)  and  (3  have  spherically  symmetric 
distributions.  As  the  lengths  of  the  left  sides  of  (2)  and  (3)  have  densities  p 
and  q  respectively,  and  the  length  of  the  left  side  of  (1)  is  r,  it  seems  that 
what  is  required  is  a  formula  for  the  density  of  the  length  of  the  convolution  of 
two  spherically  symmetric  (about  0)  distributions  given  the  densities  of  their 
lengths.  However,  as  Bickel  pointed  out  in  [1],  this  is  not  quite  correct  since 
the  use  of  the  convolution  assumes  the  independence  of  the  densities  p  and 
q,  and  these  two  densities  are  not  generally  independent,  since  in  the  right 
sides  of  (2)  and  (3)  the  terms  —  and  £m+i  are  only  uncorrelated, 
and  not  in  general  independent  unless  Eu  has  a  Gaussian  distribution. 

However,  for  large  K,  we  follow  Bickel’s  argument  in  [1]  to  show  that  the 
terms  e^f+i,!  “  ^m+i  and  £m+i  are  independent.  Let 


and 

=  (A'  —  1)  2  ^  e’Af+ij) 
i=2 

(4) 

Uk  =  £m+i,i0-  - 

(5) 

Then,  Zk  and  Uk  are  independent  and,  if  K  is  large. 

,,  (AT -1)2^  ,, 

=  Uk  ^  ^  Uk 

(6) 

and 

r.i-  K\ 

(7) 

K'^eM+\  =  \^  ^  j  +  Uk Zk- 

Thus,  for  large  K,  to  a  first  approximation  it  is  possible  to  ignore  the  depen¬ 
dence  between  (6)  and  (7).  The  terms  Uk  and  Zk  on  the  right  sides  of  (6) 
and  (7)  respectively  are  the  first  terms  in  an  asymptotic  expansion  for  large 
K  of  the  left  sides  of  those  equations. 


3.  Convolution  and  Fourier  Transform.  Let  U,  V  be  random  vectors 
in  R'^  which  are  independent  and  have  isotropic  distributions  with  densities 
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/i)  /2  respectively.  We  let  gi,g2  be  the  corresponding  densities  on  R+  of  the 
lengths  of  U  and  V,  |  [/  |,  |  V  |  .  Further,  let  /  be  the  density  oil!  +  V  and 
g  be  the  density  oi\U  +  V  \.  Then  by  the  independence  of  U,  V, 

fiy)  =  h{y  -  x)fi{x)dx.  (8) 

Our  interest  is  to  determine  a  formula  for  g  in  terms  of  gi  and  ^2-  To 
this  end,  we  introduce  the  Fourier  transform.  If  /  is  absolutely  integrable  on 
R-^,  then  the  Fourier  transform  of  /  is  defined  by 

where  the  N  vector  t  =  (ti, t^)  and  t  ■  y  =  tiyi  +  •  •  •  t^yN-  If  fiy) 
is  spherically  symmetric,  i.e.,  f{y)  is  a  function  of  r  =|  y  |  only,  say  fiy)  = 
/i(r),  then  its  Fourier  transform  fit)  is  also  spherically  symmetric;  more 
specifically,  we  have 

fit)  =  [  r^/2/i{r) J(Ar-2)/2(pr)dr,  (10) 

where  p  =|  1 1  and  Jj/(r)  denotes  the  Bessel  function  of  the  first  kind  of  order 
u.  For  a  proof  of  this  formula,  we  refer  to  Schwartz[2].  (Note  that  we  have 
introduced  the  same  letter  r  for  |  y  |  eis  we  used  for  the  density  of  the  length 
of  the  left  side  of  (1).  From  the  context,  there  should  be  no  confusion  as  to 
which  meaning  for  r  is  intended.)  The  Fourier  transform  of  the  convolution 
in  (8)  yields 

f{t)  =  hit)  flit),  (11) 

and  combining  (10)  and  (11)  gives 

/*O0  roo 

f(t)  =  P^“^y^  'rt^'^h2ir)jQ^-2)/2ipr)dr  r^^^/ii(r) J(iv-2)/2(P’")<f^,  (12) 

where  the  functions  hi,  h2  in  (12)  are  defined  by 

hiir)  =  /i(y),  h2(r)  =  /2(y),  (13) 

with  r  =1  y  I  since  the  distributions  defined  by  fx^h  are  isotropic. 

Now  if  U  is  an  isotropically  distributed  N  vector  with  density  /  on  R^, 
and  g  is  the  corresponding  density  on  R^  of  the  length  \U  \,  where  \xf'  — 
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xf,x  =  then  the  relationship  between  the  densities  /  and 

g  is  given  by 

/W  =  I  y  y  I).  (14) 

where  c^,  the  surface  ‘area’  of  the  unit,  sphere  in  R^,  is 

CN  -  .  (15) 


Hence,  upon  substituting 

hi{r)  =  fi(y)  = 

h2{r)  =  fiiy)  =  r-^^~'^'>c]^^g2{r), 


(16) 


into  (12),  we  have 

^  roo  /*oo 

f(t)=  r‘"''/^S2(r)J(K-2)/2(pr)*jf  r'^''^^gi{r)J(N-t)/2{pr)dr. 

(17) 

If  the  dimension  N  of  the  space  the  measurement  vectors  Xij  belong  to  is 
even,  say  N  =  2m,  then  J(Ar_2)/2(p^)  =  Jm-iipr)  is  a  Bessel  function  of  the 
first  kind  of  integer  order.  If  N  is  odd,  say  N  =  2m  + 1,  then  J^n-2)/2{p^)  = 
Jm-i/2{p^)  is  a  Bessel  function  of  the  first  kind  of  firactional  order,  and  is 
closely  related  to  the  Spherical  Bessel  function  of  the  first  kind  jn{z),  defined 
by 

jn{z)  =  lJ^Jn+l/2{z).  (18) 

More  detailed  information  about  Bessel  functions  may  be  in  found  in  Abramowitz 
and  Stegun  [3]. 


4.  Convergence  of  the  Integrals,  We  now  discuss  convergence  of 
the  integrals  in  (17).  Both  integrals  are  functions  of  the  variable  p  =|  t  |, 
hence  f{t)  is  a  spherically  symmetric  function  of  the  transform  variable  t. 
This  means  that  the  inverse  Fourier  transform  of  f{t),  i.e.  /(y),  can  also  be 
obtained  as  a  one-dimensional  integral  with  a  Bessel  function  kernel. 

For  fixed  p  and  small  r,  J(n-2)/2{p^)  =  0(r(^"'^)/^).  Since  we  expect 
gi^\o)  =  0  for  z  =  1, 2  for  0  <  j  <  iV  —  1  (see  Bickel  [1]),  it  follows  that  both 
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integrals  in  (17)  are  convergent  at  the  lower  limit  0.  At  the  upper  limit  oo, 
J{N-2)/2{p^)  =  and  this  by  itself  will  not  be  enough  to  make  the 

integral  convergent.  However  the  term  will  also  make  the  integrals 

converge  at  the  upper  limit  if  N  is  sufFciently  large.  In  practice  the  densities 
9x^92  also  tend  to  zero  sufficently  rapidly  to  make  the  integrals  converge  at 
the  upper  limit.  If  these  densities  have  compact  support  (i.e.,  are  zero  outside 
of  a  closed  bounded  subset  of  R"*"),  then  the  integrals  in  (17)  no  longer  have 
infinite  upper  limits.  The  integrals  in  (17)  can  be  evaluated  accurately  and 
efficiently  by  standard  numerical  quadrature  methods. 


5.  Inversion.  Note  that  the  right  side  of  (17J)  is  a  function  of  p  =|  t  | 
only.  Hence  f{t)  is  spherically  symmetric.  Let  f{t)  —  G{p).  The  inverse 
Fourier  transform  of  f{t),  i.e.,  f{y)  from  (9),  is  defined  by 

It  follows  by  analogy  with  (10)  that 

=  f  p^^^G{p)J(^N-2)/2{pr)dp,  (20) 

JO 

where  r  =|  y  |.  The  relationship  (14)  between  the  density  /  of  an  isotropi¬ 
cally  distributed  N  vector  and  the  corresponding  density  g  of  its  length  then 
implies 

roo 

g{r)  =  I  p^f'^G{p)J(N-2)/2{pr)dp  (21) 

where  G{p)  =  f{t)  is  given  by  (17).  Numerical  evaluation  of  the  integral 
in  (21)  proceeds  similarly  to  the  integrals  in  (17).  With  gx,g2  taken  as  p,q 
introduced  at  the  end  of  section  1,  and  g  taken  as  r  (the  density  function  for 
I  Xm+i,x  ~  I)  woi  I  y  I),  (17)  and  (21)  together  give  the  density  r  in  terms 
of  p  and  q. 


6.  Conclusions.  We  have  shown  that  for  large  K,  the  left  sides  of 
(2)  and  (3)  are  approximately  independent.  Under  this  approximation,  we 
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derived  an  expression  for  the  density  r  of  the  length  of  the  left  side  of  (1)  in 
terms  of  the  densities  p,  q  of  the  lengths  of  the  left  sides  of  (2)  and  (3).  This 
result  is  contained  in  equations  (17)  and  (21)  of  the  previous  section. 

Future  work  could  include  determining  corrections  to  the  above  result  for 
a  finite  number  of  measurements  K,  and  practical  numerical  implementation 
of  the  above  result. 
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